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Abstract 

Possible mixing of the active and sterile neutrinos has been considered both in the 
single Am 2 approximation and in the case of more than one Am 2 . We perform global 
fit of the available solar neutrino data with free boron neutrino flux in the single Am 2 
context. The best fit value corresponds to zero fraction of sterile component g = 0. We 
get the upper bounds: g < 0.26(0.64) at la (3a). Due to degeneracy of parameters no 
one individual experiment restricts g. The bound appears as an interplay of the SNO 
and Gallium as well as SuperKamiokande data. Future measurements of the NC/CC 
ratio at SNO can strengthen the bound down to < 0.5 (3c). If KamLAND confirms 
the LMA solution with its best fit point a combined analysis of the KamLAND and 
solar neutrino results will lead to g < 0.19(0.56) at la (3a). We find that existence of 
sterile neutrino can explain the intermediate value of suppression of the KamLAND 
event rate: Rkl ~ 0.75 — 0.90 in the case when more than one Am 2 is involved. 


1 Introduction 


Sterile neutrino, u s , is a panacea of various problems which appear from time to time in 
neutrino physics. Oscillation interpretation of the LSND result [1] in the context of schemes 
which describe the solar and atmospheric neutrino problems [2] is the long term motivation 
for existence of v s . 

There are, however, other motivations to search for the sterile component. It might 
happen that v s is the right-handed component of neutrino which turns out to be light (or 
massless) in the context of the see-saw mechanism, as a consequence of certain symmetry [3]. 
Such a symmetry (L e — L fJ — L T as in [3] ) can also be responsible for large or maximal lepton 
mixing. Sterile neutrino can be a component of multiplet of the extended gauge group, e.g., 
Eq [4]. Sterile neutrino - light singlet fermion - can originate from completely different sector 
of theory. Many extensions of the standard model predict existence of singlet fermions: u s 
can be a mirror neutrino [5], goldstino in SUSY [6], modulino of the superstring theories [7], 
bulk fermions related to existence of extra dimensions [8], etc.. 

Neutrinos are unique particles, since only they can mix with singlets of the Standard 
Model. This mixing can lead to coherent effects which provide high sensitivity to existence 
of v s . The mixing terms as small as (10~ 10 — 10 -6 ) eV can produce the observable effects 
and in the case of supernova neutrinos the masses can be even smaller. In a sense, neutrino 
mixing is a window to physics beyond the SM. 

Another motivation to search for and restrict sterile component is that interpretation 
of certain results can be changed strongly if one admits even very small admixture of sterile 
neutrinos. This concerns with the neutrinoless double beta decay, CP-violation, generation 
of maximal/large flavor mixing [9], etc.. 

So, we deal with long term program of searches for sterile neutrinos and improvements 
of the bounds on mixing of sterile neutrinos with all possible masses. 

In this paper we will consider in details searches for sterile component in the solar 
neutrino flux. First detailed discussion of signatures of the sterile neutrinos have been done 
in [10, 11, 12] where it was marked, in particular, that studies of the spectrum distortion 
can reveal existence of v s . 

There were a number of fits of the solar neutrino data in terms of pure v e —» v s sterile 
conversion, but recent SNO results [13, 14] exclude pure v e —» v s conversion at high con¬ 
fidence level. Still partial transformation of v e in v s is possible. The effects, and allowed 
fraction of sterility depend on specific scheme of mixing. Most of recent studies have been 
done in a single Am 2 context: according to which the electron neutrino mixes with the state 
being the combination of the active (u /t and u T ) and sterile neutrino. Fraction of sterility 
is described by the parameter rj. In [15] it was marked that p is weakly restricted provided 
that the original boron neutrino flux is substantially larger than the SSM flux. Global fit of 
the solar neutrino data [16] gives rj < 0.25 (0.6) at la (3a) level. If the boron neutrino flux 
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is fixed according to the SSM with corresponding uncertainties the bound becomes stronger 
at the 3a level: rj < 0.5 [17]. 

In this paper we will consider effects of sterile neutrinos and bounds on sterility from 
the existing and future experiments. In particular, we study how combined analysis of the 
KamLAND and solar neutrino data can improve the bound. We generalize the analysis for 
the case when more that one Am 2 is relevant. We discuss how presence of sterile component 
can change predictions for KamLAND rate. 

The paper is organized as follows. In section 2 we consider the mixing the sterile neutrino 
in a single Am 2 context. In section 3 we perform global fit of all solar neutrino data, and 
put bounds on the sterile fraction rj. In section 3.2 we identify the data sensitive to rj. In 
section 4 we study how future experiments can improve the bound. We consider precision 
NC measurements at SNO (section 4.1) and combined analysis of KamLAND and solar 
neutrino data (section 4.2). In section 5 we consider the v s mixing in the case when more 
than one Am 2 is relevant. We discuss how the presence of sterile component can change 
predictions for the KamLAND experiment. We present our conclusions in section 6. 


2 Mixing of sterile neutrinos 

In what follows we will consider the case of one sterile neutrino. In general, fourth neutrino 
may have an arbitrary mass and sterile component can mix with all three active neutrinos. 
So, the scheme will have 4 new real parameters: the mass m 4 and 3 mixing angles. Situation 
can, however, be simplified if one takes into account existing bounds on mixing of active 
neutrinos. We will consider the 4-neutrino schemes which explain the solar and atmospheric 
neutrino data. We assume that two mass eigenstates, v\ and u 2 , are splited by the solar 
Am 2 < several x 10 -4 eV 2 . We will call them the “solar pair”. Then the state u 3 , is splited 
from the solar pair or from /z 4 by the atmospheric mass split: Am 2 atm ~ (2 — 3) • 10 -3 eV 2 . 
We define the mixing matrix, 11 U m \ | as the unitary matrix which connects the flavor and 
the mass eigenstates: v a = J2i ce = e, i = 1, 2, 3,4. 

2.1 Single Am 2 case 

As far as solar neutrinos are concerned, mixing of the electron neutrinos plays crucial role. 

Let us first assume that U e3 = U ei = 0, so that the electron flavor is distributed only 
in the “solar pair of states”: In fact, there is a strong upper bound on U e3 from CHOOZ 
experiment [19], and strong bound on I / e 4 from BUGEY experiment [20] if corresponding 
Am 2 is large. 

Since only U e \ and U e2 are non-zero the orthogonality conditions for v e and other neu- 
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trinos, v tll v T , v s can be written as 


uhu al + U} 2 U a2 = 0 , a — s, h,t . 


From these equations we get immediately: 

Uj* = Ut 2 = Us2 = _uU , . 

u hi u rl u s i uh' 1 J 

That is, all non-electron flavor components enter zq and zq in the same combination. Indeed, 
let zq be the combination of the non-electron neutrinos which mixes with v e in zq : 

1 

Vl = cos d X2 v e - sin e 12 u x , zq = V U al v a , (3) 

Jl-U^ a 


where sin Q\ 2 — y l — . Similarly, let v y be the combination of the non-electron neutrinos 

which mixes with v e in zq: 


zq = sin d 12 zq + cos # 12 ^, zq = 


'l - I/. 2 , 


'y ) U a2 V a . 


Then according to (2), we have: 


v x = v y . 


Thus, the scheme is reduced exactly to the two two neutrino case: v e mixes with zq in the 
states v\ and v 2 and the mixing angle equals $ 12 - 

The state v x can be written as 

v x = y/1 — ‘Wa + (6) 

where v a is the combination of active (non-electron) neutrinos, zq and zq: 

1 

Va / —(h^/1 V y T U T \V T ^j. (7) 

7 1 - 

So, y/rj describes admixture of sterile neutrino in the state v x to which v e can be transformed. 
According to (3), (4) and (6), admixtures of the sterile component in zq and v 2 equal 


Usi = ^rjsin 9i2, U s2 = ^/rjcos d 12 , 


and consequently, 

\U s1 \ 2 + \U s2 \ 2 = V- (9) 

That is, r] gives total fraction of the sterile component in the solar pair. 

We have proven that if the electron neutrino is distributed in two mass eigenstates 
only, then independently of other mixings and masses, the effect of sterile component in 
solar neutrinos is described by one parameter only, which is the total amount of sterile 
component in these two mass eigenstates. 
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2.2 Conversion probabilities and observable fluxes 


Let us consider the oscillation effects. We introduce 


Pee = Pee(Am 2 u , 0 12 , V ex ) (10) 

- the v e survival probability in the system of mixed neutrinos v e —v x . Here V ex = V e —V^{l—rj) 
is the effective potential. According to (3) and (6) the transition probabilities of the electron 
neutrino to the sterile, P es , and active, P ea , components in terms of P ee equal 

P es = 77(1 - P ee ) 

P ea = (1 — 77)(1 -P ee ). (11) 


Using these probabilities we can write the fluxes of neutrinos which determine the rates of 
events of different types. Let us introduce fa, the original boron neutrino flux in units of 
the SSM flux, : 


Ib = 


B 


pSSM 


( 12 ) 


Here the SSM boron neutrino flux is taken to be F^ Shl = 5.05 • 10 6 cm -2 s~L Then the 
clastic scattering events (ES) (detected by SuperKamiokande and SNO), the neutral current 
(NC) and charged current (CC) event rates at SNO are determined by the following fluxes: 


^cc — fsPee, (13) 

$NC = f B \Pee+(l-V)(l-Pee)}, (14) 

$ES = f B \Pee + r(l-ri)(l- Pee)}, (15) 


where r = 0.16 is the ratio of the i/ M)T e and z/ e e elastic scattering cross sections. The 
equations (14) depend on three parameters, f B , Pee and r) via two combinations 


x = f B Pee, y = f B (l-rj)(l-P ee ). (16) 

In terms of variables x and y , the fluxes can be rewritten as 

= x + ry, 4>cc = x, $ NC = x + y. (17) 

Excluding x and y we get relation between the fluxes [15]: 

®es = 4>cc + 'KdWc - ®cc) = *hcc( 1 - r) + r<h NC (18) 

which has a simple interpretation (see the first equality): the flux measured in the neutrino- 
electron scattering equals the flux of electron neutrinos at the detector plus the flux of non- 
electron active neutrinos, &nc ~ $cc, suppressed by r. The equality (18) is well satisfied 
in the experiment. 
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The equalities (17) imply the degeneracy of parameters [15]: different values of the 
original parameters /#, 77 , P ee lead to the same observables provided that they keep to 
be constant the combinations x and y. In particular, changes of r) can be compensated 
by corresponding variations of P ee and fs, so that the combinations and therefore the 
observables are not changed. Clearly this is possible if P ee does not depend on the energy or 
time. Therefore the energy and the time variations of the probability break the degeneracy. 

In other words, time variations and distortion of the spectrum are, in general, sensitive to 
the fraction of sterile neutrinos. Since no significant distortion or variations have been found 
experimentally, the present data have no strong sensitivity to the fraction 77 as we will see 
from exact calculations in the sect. 3. 

3 Global fit and bounds on fraction of sterile compo¬ 
nent 

3.1 Global analysis 

We have performed a global analysis of all available solar neutrino data taking into account 
possible presence of the sterile component. We follow the procedure described in our previous 
publication [ 21 ] where details can be found, and here we summarize the main ingredients of 
the analysis. 

We use the same data sample as in [21], which consists of 

- three total rates: 1) the Ar-production rate, QAr, from Homestake [22], 2) the Ge— production 
rate, Qc e from SAGE [24] and 3) the combined Ge—production rate from GALLEX and 
GNO [25]; 

- 44 data points from the zenith-spectra measured by Super-Kamiokande [26] during 1496 
days of operation; 

- 38 day-night spectral points from SNO [13, 14]. 

All solar neutrino fluxes (but the boron neutrino flux) are taken according to SSM 
BP2000 [28]. We use the boron neutrino flux as free parameter. For the hep—neutrino flux 
we take fixed value Fh, ep = 9.3 x 10 3 cm -2 s _1 [28, 29] . 

In contrast to [21], the neutrino scheme includes now mixing with sterile component 
which is described by the parameter 77 ( 6 ). So, there are three oscillation parameters: the 
mass squared difference, Am 2 , and two mixings: tan 2 6 and 77 . Consequently, in the free 
boron neutrino flux fit we have four parameters: Am 2 , tan 2 9, 77 and /#, and therefore there 
are 81 (data points) - 4 = 77 d.o.f.. 

We perform the y 2 test of the oscillation solution by calculating 

^global = Xrate + XsK + XsNOi ( 19 ) 
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V 

0.0 

0.2 

0.4 

0.6 

x 2 

65.1 

65.5 

67.5 

72.6 

Am 2 (xlO -5 ) 

6.3 

5.5 

4.6 

4.0 

tan 2 9 

0.41 

0.35 

0.30 

0.26 

,1b 

1.046 

1.188 

1.356 

1.570 


Table 1: Minimum y 2 for different values of n. Also presented are the values of parameters 
Am 2 , tan 2 9 and rj that minimize y 2 - 


where y 2 ate , Xsk an4 Xsno are ^ ie contributions from the total rates, the Super-Kamiokande 
zenith spectra and the SNO day and night spectra correspondingly. Each of the entries in 
(19) is the function of four parameters (Am 2 , tan 2 6, f B , rj), in particular, 

X 2 global = X 2 global( A , tan 2 9, f B , fj) . (20) 

We perform the global analysis for various fixed values of rj. We first find the best fit 
points for a given rj minimizing x^iobai with respect to Am 2 , ta n 2 9 and f B . This gives 
Xminiv)- The results are shown in the Table 1. The best fit corresponds to zero value 
of sterile fraction and f B = 1.05. So, the absolute minimum is characterized by Xmm(0)- 
According to the Table 1, x 2 increases with rj, and moreover, the increase is fast for rj > 0.3. 
Ay 2 = Xmm( 0 - 4 ) - Xmm(°) = 2 - 4 ; for V = °- 6 this difference is Ax 2 = 7.5. 

Then we construct the 3a regions in the Am 2 — tan 2 9 plane (Fig. 1) minimizing 
X 2 (A m 2 , tan 2 9, f B ) with respect to f B : we denote the corresponding minimal value of x 2 by 
X 2 (Am 2 , tan 2 6). Then 3a contours are determined by y 2 (Am 2 , tan 2 9) — Xminiv ) = H-73, 
where Xminiv) i s the minimum x 2 for a given value of rj. We have found also lines of constant 
values of f B which minimize x 2 f° r a given rj. 

According to Fig. 1 and the Table 1, with increase of rj the allowed region shifts to 
smaller values of tan 2 9 and Am 2 . At the same time f B and y 2 nm increase. 

These features can be well understood using Eqs. (15 - 16). Increase of rj [decrease of 
(1 — rj)] can be compensated by increase of f B and decrease of P ee . Since P ee ~ sin 2 9, the 
decrease of P ee implies decrease of mixing angle. With decrease of mixing the distortion 
of the spectrum (mainly the turn up of the probability at low energies) becomes stronger. 
This leads to shift of the allowed region to smaller values of Am 2 (shift of the spectrum 
from the adiabatic edge). 

Let us consider breaking of parameter degeneracy. Taking P ee « sin 2 9 in expressions 
(16) and using results of the Table 1 we find values of parameters x{rj) and 2 /( 77 ) for different 
rj in the best fit points (see Table 3.1). As follows from the Table 3.1 in the best fit points 
the parameters x{jj) and 7 /( 77 ) are not invariant under changes of rj: x slightly increases (by 
5%) whereas y decreases significantly: almost by 50%. Another way to see these variations 


6 



V 

0.0 

0.2 

0.4 

0.6 

X 

0.304 

0.308 

0.313 

0.320 

y 

0.742 

0.703 

0.623 

0.492 

y(x — const) 

0.726 

0.694 

0.622 

0.500 


Table 2: Values of parameters x and y in the best fit points of global analysis for different 
values of rj. 


is to keep x = $cc = constant and to exclude P ee from the combination y: 

y=( 1 ~v)(fB-$cc)- (21) 

Using $>cc — 0-32 we find y(rj) shown in the third line of the Table 3.1. Thus, the invariance 
of x and y (and therefore a degeneracy) is broken which leads the dependence of Xmin on V 
and therefore to sensitivity of the analysis to the presence of sterile neutrinos. 

For fixed values of rj and Jb we minimize y 2 with respect to Am 2 , tan 2 6 which gives 
Xa 0 = XAe( ,r h fs)- Using this function we construct contours of constant confidence level 
in rj, fs plane (fig. 2). From this figure we get an upper limit for the sterile fraction 

rj < 0.38 (0.70), la (3 a) (22) 

The boron neutrino flux which correspond to maximal allowed value of rj equals fs = 1-25 
(la) and fs = 1-68 (3a). 

We have performed minimization of y 2 with respect to Am 2 , tan 2 6 and fs, thus finding 
y a,6»,/ b = X 2 (v)- Then the function y 2 ( 77 ) allows to put the bounds for 1 d.o.f.: 

y < 0.26 (0.64), la (3a) (23) 

These results are very similar to results obtained in [16, 17]. In particular, it was found [16] 
that at la: rj max = 0.25 for fs ~ 1.25, and at 3a: rj max = 0.6 for fs ~ 1.75. 

Imposing the SSM bound on the neutrino flux further strengthen the bound on “steril¬ 
ity”. It was found in [17]: rj < 0.25 at the la and rj < 0.5 at the 3a level. Notice that in 
fact la bound is unchanged since at this level the required boron neutrino flux is within la 
SSM prediction. In contrast, the limit becomes stronger at the 3a level, where substantially 
larger than in SSM flux is needed. 

For the SSM value of the boron neutrino flux, fs = 1, we get from the fig. 2 the following 
bounds on the sterile fraction: rj < 0.14 (la), rj < 0.30 (2a), rj < 0.47 (3a). 

3.2 Who does not like sterile neutrino? 

Let us identify observables which are sensitive to rj. As we have discussed in sect. 2.1 
these observables should contain information about the time variations or/and the energy 
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dependence of the conversion effect. These dependences remove the degeneracy and therefore 
restrict rj. 

No statistically significant time variations have been found although some indications 
of the Day- Night asymmetry, A ND , exist [26, 14]. Let us consider if A ND can restrict 
the admixture of sterile neutrino. The high energy neutrino data can be described by the 
zvsurvival probability during the day and during the night: P^, PjJ, or equivalently, by 
A P ee = P^ e — P® and P ee = (P e ^ + P®)/2. Consequently, in Eq. (13 - 15) we should 
substitute P ee P ee and add to analysis an additional observable: 

And = —jr 1 - ( 24 ) 

* ee 

In the LMA region the average suppression is determined basically by mixing angle: P ee ~ 
sin 2 6 , whereas the DN asymmetry depends mainly on Am 2 . This can be seen from lines of 
constant ratio CC/NC ~ P ee and lines of constant And which are nearly orthogonal each 
other (see figs. 5 and 6 from [21]). Therefore A P ee and P ee do not correlate and can be 
considered as independent fit parameters. For a given P ee , the difference A P ee can change in 
the wide range, and the data on the asymmetry can be easily fitted. So, adding experimental 
information about the Day and Night asymmetry does not help to resolve degeneracy and 
therefore to improve the bound on rj. Precise measurements of the zenith angle dependence 
of signal will change this conclusion. 

No statistically significant distortion of the energy spectrum is found for E > 5 MeV, 
which supports the degeneracy problem. The distortion (in experiments sensitive to E > 5 
MeV) is expected for higher values of Am 2 and small mixing. However, substantial distor¬ 
tion is expected over wider detectable energy range which includes also the pp— neutrinos. 
Therefore, the interplay of the high energy data and results from the low energy experiments 
should break the degeneracy. 

To illustrate this effect, let us describe the distortion by two (for simplicity) different 
values of probabilities: at high energies, Pand at low energies, P^ e (.E < 0.5 MeV). 
The signal in the Gallium experiment is determined essentially by P^ e : so that the (Re¬ 
production rate equals Qa e °c P e ^, and obviously, no substantial dependence on fs appears. 
In contrast to the Day-Night asymmetry case, the probabilities P^ e and P^ e are strongly 
(anti) correlated. Indeed, in the range of pp-neutrinos: 

P e ^~ 1-0.5 sin 2 26, (25) 

and since P^ e ~ sin 2 6 we get 

P e L e ~l-2P e H e (l-P e H e ). (26) 

With decrease of P the low energy probability P^ e increases. Consequently, the predicted 
value of the Ge-production rate increases. Thus, combination of the Gallium and SNO data 
should break the degeneracy of parameters. 



V 

0.0 

0.2 

0.4 

0.6 

x 2 

25.6 

25.4 

25.3 

25.1 

Am 2 (xlO -5 ) 

4.5 

4.5 

4.0 

3.6 

tan 2 9 

0.45 

0.35 

0.28 

0.18 

1 b 

1.039 

1.201 

1.432 

2.013 


Table 3: Minimum Xsno f° r different values of rj. Also presented are values of parameters 
Am 2 , tan 2 9 and rj that minimize the Xsno- 


Now let us describe results of exact numerical calculations. 

1) . We have performed the x 2 analysis of the SNO data only. The best fit points in the 
oscillation plane as well as /# which minimize X 2 sno( 7 i) are given for different values of rj 
in the Table 3.2. Notice that with increase of rj the quality of the fit does not change: y 2 
even slightly decreases. So, as far as SNO data alone are concerned, no bound on rj appears. 
Variations of rj can be completely compensated by changes of P ee and /#. The day-night 
asymmetry has low statistical significance. 

The SNO experiment has higher sensitivity to the NC- events. For this reason the 
allowed region provided by SNO is restricted in tan 2 9 from both sides. Increase of rj has 
the effect of moving the allowed region to smaller values of tan 2 9, as can be seen in fig. 3. 

Indeed, to keep the combinations (16) unchanged with increases of rj one needs (i) to 
diminish P ee , which implies decrease of 6 since in the LMA region P ee ~ sin 2 9, (ii) to 
decrease Am 2 to avoid distortion (turn up) of the spectrum at small energies and (iii) to 
increase J'b to compensate decrease of P ee . 

2) . In fig. 3 we show also lines of constant Ge-production rate. The rate increases with 
decrease of mixing. According to the figure, in the SNO allowed region for rj = 0.6 the rate 
Qce > 75 SNU and in the best fit point Qc e — 83 SNU. For rj = 0.4 we get Qc e = 75 SNU 
in the best fit point. The combined results from SAGE, Gallex and GNO experiments is 
Qce = (70 ±4) SNU. So, the Gallium results prevent from further shift of the allowed region 
to smaller 9 12 and therefore forbid larger rj. To illustrate the role of Gallium experiments, 
we find the bounds on rj that can be derived by taking only SNO and Gallium data: 

rj < 0.43 (0.78), la (3a), (27) 

which are weaker than those obtained by taking all the data (23). 

In principle, further more precise measurements of the Ge-production rate could im¬ 
prove the bound on the sterile component. 

3) . Let us comment on the impact of SuperKamiokande. In fig. 4 we show result of 
analysis of the SuperKamiokande data only for different values of rj. The excluded (at 
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la) regions in the oscillation parameter space are clue to the effects of spectrum distortion 
and the Earth v e regeneration. They are only slightly modified by the presence of sterile 
component. With increase of 77 the excluded region shifts to larger mixing angles. This 
is related to the fact that with increase of sterile fraction the distortion of the spectrum 
increases due to decrease of the damping effect related to active neutrinos. The exclusion 
due to earth regeneration becomes weaker. Since for active-sterile system the effective 
matter potential is smaller the regeneration region shifts to smaller Am 2 . Correspondingly, 
lower bound of the allowed region shifts to smaller Am 2 . 

In the allowed region, where distortion and time variations are small the data can be 
reproduced for different values of r / by adjusting fs and P ee (see (15)). 

In fig. 4 we show also fs that gives the best fit to the Super-Kamiokande data. With 
increase of 77 the lines fs = const shift to larger tan 2 6 and Am 2 (where the survival 
probability is larger) to compensate decrease of muon and tail neutrino contribution. Large 
values of rj are allowed, when the boron neutrino flux is much larger than the SSM flux. 

Similarly to the Gallium experiments, SuperKamiokande prevents from the shift of the 
allowed region to smaller mixings. Consequently, the combined fit of the SNO and SK will 
also give the bound on the sterile fraction. Comparison of the bounds (23) with and (27) 
without the SK result allows to evaluate the role of SuperKamiokande. 

4). The Homestake experiment alone is not sensitive to rj. However, it has some im¬ 
pact since its result restrict the oscillation parameters. The Hr-production rate decreases 
with mixing angle and the agreement with the Homestake result becomes even better with 
increase of 77 . Comparing the SNO and Homestake results we find that in the best fit SNO 
point for 77 = 0.6 the rate equals Qa t = 2.75 SNLI which is within la of the Homestake 
result. So, the SNO and Homestake data favor presence of sterile component. 


Summarizing, the bound on sterility appears as combined effect of several experiments. 
In the global analysis the constrain on 77 comes from the SNO measurements of the spectrum 
which contains information about the NC and CC event rate from the one hand side and 
Gallium as well as the SK experiment from the other side. With increase of 77, SNO pushes 
the allowed region to smaller tan 2 6, whereas Gallium and SK experiments put lower bound 
on tan 2 d. For instance, the best SNO fit point for 77 = 0.6 lies in the excluded Gallium 
region and at the border of the SK excluded region. 
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4 


Future experiments and bounds on 77 


4.1 NC measurements at SNO 

SNO is taking data with much higher sensitivity to the NC current events, and therefore 
precision of measurements of &nc will be significantly improved. Clearly, distortion of the 
spectrum of the NC events will be signature of the sterile neutrino conversion. Although if 
LMA is the dominant solution, one does not expect observable effect. As far as total rate is 
concerned, the parameters can be adjusted in such a way that even more precise data can 
be reproduced for large rj. Basically we will get the same shift of the allowed region as in 
fig. 1 but the size of the region around the the best fit point will be smaller. For this reason 
no drastic improvements of the bound on sterility is expected. 

To evaluate the impact of future WC-measurements at SNO we show in fig. 5 the lcr 
region of predictions of &nc as a function of rj from the global fit. The results have been 
obtained in the following way. We use the global analysis of the solar neutrino data for 
different values of rj, as is shown in fig. 1. Then we calculate &nc in the best fit points 
(solid line) as well as the spread of predictions for &nc within la allowed region. Using 
the fig. 5, we find that the present la experimental result on ®nc, gives the upper bound 
rj <0.7 which agrees with the 3<r bound from the global analysis. Assuming three times 
smaller error bars for the NC measurements, <f> 7 vc = 1 ± 0.03, we get the bound rj < 0.5 
which should be compared with present bound rj < 0.64. 

4.2 KamLAND and sterile fraction 

The signal at KamLAND is determined by the vacuum oscillation survival probability, so 
that the rate of events is 

2 FiaA • < 28 > 

where F % is the flux from i reactor, a is the cross-section of up —» e + n reaction, A is the 
energy resolution function. The sum is taken over all reactors contributing to the flux at 
Kamioka and the integral denotes schematically the integrations over the neutrino energy, 
and the true energy of produced positron. This rate does not depend on rj or /#. Being 
combined with the solar neutrino results it breaks degeneracy. KamLAND will fix the 
oscillation parameters A m\ 2 , tan 2 # 12 , and consequently, will allow to predict P ee which 
describes the solar neutrino conversion. 

Let us evaluate how combined analysis of the solar data and KamLAND can improve 
the bound on sterile component. 




dE 1 - sin 2 26 
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Let us introduce the suppression factor of total event rate above certain threshold: 


R KL (Am 2 , tan 2 9) 


N(Am 2 , tan 2 9) 
iVo 


( 29 ) 


where N is given in (28) and N 0 is the event rate in absence of oscillations. 

Since predictions for KamLAND do not depend on the state to which the electron 
neutrino converts, we can use results of calculations of R KL for the pure active case [21]. In 
fig. 8 of paper [21] we have presented lines of constant Rkl hr (Am 2 , tan 2 9) plane. In the 
best fit points we find Rkl( 0) = 0 . 66 , Rkl( 0 - 2 ) = 0.63, Rrl( 0-4) = 0.55, Rkl( 0 - 6 ) = 0.52. 
Due to decrease of Am 2 in the best fit points, the expected value of the suppression in the 
KamLAND experiment becomes stronger. 

To quantify impact of the event rate measurements at KamLAND we have constructed 
the \ 2 (Am 2 , tan 2 9, f B , rj) as the function Rkl- Practically we impose the following condi¬ 
tion Rkl( Am 2 , tan 2 9) = const which implies relation between Am 2 and tan 2 9. So, for a 
given value Rkl the y 2 is the function of Am 2 or tan 2 9 only: 


x 2 = X 2 (tan 2 6, Rkl, !b, rj). 


(30) 


For each value of Rkl we minimize y 2 with respects to tan 2 9, f B , rj thus finding x^Js^Rkl)- 
This function is shown in the upper panel of fig. 6 . In the bottom panel the corresponding 
values of f B and r/, which minimize y 2 , are presented. For comparison in the upper panel 
we show also the function X0j b (Rkl) for V = 0. 

The following remarks are in order. 

For R K l < 0.67 the quality of the fit does not depend on Rkl- It corresponds to 77 = 0 
and fs ~ 1.05. For Rkl > 0.67, y 2 increases with Rkl- Simultaneously the best fit values 
of r] and fs increase. At 3 a level (y 2 = 74) we find maximal value of Rkl and corresponding 
values of 77 and f B 

Rkl ~ 0.75, 77 = 0.35, f B = 1.4. (31) 

For 77 = 0 the value Rkl ~ 0.725 is allowed. So, the introduction of sterile neutrino in the 
single Am 2 approach does not lead to significant increase of maximal allowed value of Rkl- 
The goodness of the fit decreases sharply: e.g. for Rkl = 0.8 we get Ay 2 = 25 even in 
presence of sterile component. 

Let us analyze how the KamLAND spectral measurements can help in constraining 
77 . We have simulated the KamLAND spectral data for 4 different points a), b), c), d) 
in the Am 2 — tan 2 9 plane (marked by stars in fig. 7a)). In simulations we use statistics 
which can be collected during 3 years in 600 tons above the energy threshold for visible 
energy E vls = 2.6 MeV. In absence of oscillations this would correspond to 1170 events. We 
assumed a 4% uncertainty in the overall normalization and a 2% uncertainty in the energy 
scale. We have simulated the energy spectrum (for each point) using 12 bins of 0.5 MeV 
size in the visible energy. 
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We perform y 2 analysis of simulated spectra assuming a full correlation of systematics 
uncertainties. We minimize 

Xkl = /A Ri a i,j Rj (32) 

i=l,12 1=1,12 

and construct contours of constant la, 2a, 3a confidence level around selected points in 
the Am 2 , tan 2 9 plane (see fig. 7a). Notice that our allowed regions are somehow larger 
than those obtained in [31, 32, 35], since we include systematic errors and also use higher 
threshold (smaller number of events). In agreement with results of calculations of other 
groups [31, 32, 35] we find that KamLAND will be able to measure Am 2 rather precisely 
(if Am 2 < 2 ■ 10” 4 eV 2 ). At the same time determination of the mixing angle will not be 
very accurate. For the present best fit point (5) we get from the figure tan 2 d = 0.28 — 0.60 
at the la level, and tan 2 9 = 0.25 — 0.85 at the 3a level. 

Then we have performed combined analysis of the existing solar and simulated Kam¬ 
LAND data (notice that by the time KamLAND will have 3 years statistics new solar 
neutrino data will appear). For each selected point of the oscillation parameters in fig. 7a 
we calculate a global y 2 ot , which includes contributions from the solar neutrino analysis, 
xiw and the KamLAND, Xklv 

Xtot = Xkl ( tan2 0, Am 2 ) + y| un (tan 2 9, Am 2 , rj, f B ) ■ (33) 

For each value of rj we minimize Xtot with respect to Am 2 , tan 2 6 and f B - Results of this 
minimization as a function of rj are presented in fig. 7b. Let us comment on results for four 
selected points. 

The point b ) is the best fit point of the global fit which corresponds to rj = 0. As follows 
from the figure, 

rj < 0.19 (0.56) la (3a), (34) 

thus improving the bound (23) from the solar data analysis only. 

The points c) and d) are at larger Am 2 and tan 2 612 as compared with b) thus further 
removed from the region preferred in the case of non-zero sterile component (small Am 2 
and tan 2 $ 12 ). If KamLAND selects these points the bound on sterile component will be 
stronger than in a): rj < 0.09 (0.45) at la (3a). 

The point a) is in the region preferable by sterile component. Selecting this point 
KamLAND will favor non-zero rj. From the lower panel we get: rj = 0.36io;}2 a ^ 1°" and 
rj > 0 at 2a level. The la lower bound is rj > 0.24. 

For the points a) - d) we get the upper bounds on rj: rj < 0.08 (c and d), 0.19 (b) and 
0.52 (a). We have performed similar analysis for all the points of the oscillation parameter 
space. We show in fig. 8 the results as lines of constant upper bounds in the oscillation 
parameter plane. As a tendency, the bound becomes weaker with decrease of mixing and 
Am 2 . For sufficiently small tan 2 9 the best fit corresponds to a non-zero value of rj and the 
lower bound on rj appears. In the lower panel of fig. 8 we show such a lower limit for all 
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points in (Am 2 , tan 2 9) plane. 


5 Beyond single Am 2 context 

The results described in the previous sections can change if the electron neutrino mixes in 
the other mass eigenstates beyond the solar pair. Now more than one Am 2 is involved, and 
effect depends on specific values of additional Am 2 . 

In what follows we will assume that mixing of v e in the states z /3 and z / 4 is small. Then 
if additional Am 2 are outside the range of the MSW enhancement (e.g. Am 2 4 is large) the 
effect of new v e mixing is small being suppressed by the mixing itself. If, however, Am 2 4 
is in the range of the MSW enhancement, the effect of additional mixing can be large. We 
will consider both possibilities. 

5.1 Large additional Am 2 

Suppose v e has non-zero admixture in v 3 characterized by by U e 3 = S 13 . (Recall that z /3 is 
separated from the solar neutrino pair by the mass gap which corresponds to the atmospheric 
neutrino split.) We assume that ly-adrriixture in z / 4 is negligible so that the mass of z / 4 can 
be arbitrary within the allowed range by non-solar neutrino experiment. We assume also 
that in the limit 543 — > 0, the state z /3 coincides with v' a - the combination of ry and z/ T , 
which is orthogonal to v a (7). In turn, z / 4 is the combination of u s and v a . As a result, we 
get the following mixing pattern: 

V\ = cos 6*12 (ci 3 z/ e — Si 3 v' a ) — sin 6 * 12 ^, 
v 2 = sin 0 i 2 (ci 3 i/ e - s 13 u' a ) + cos 012 V*, 

*4 = c 13 u' a + s 13 u e , 

z / 4 = *4, (35) 

where v x is defined in ( 6 ) and v' x is the orthogonal to v x combination of v s and v a . The 
mixing (35) differs from the generic case since v e and v' a do not mix in z/ 4 , and u a and u s do 
not mix in z/ 3 . 

The dynamics of propagation is similar to that of the case of three neutrinos. Indeed, 
let us represent z/ e in the following form: 

*4 = ci 3 z/' + S 13 Z /3 , (36) 

where z/ = cos 0 i 2 Z/i +sin 6 l 1 2 Z /2 coincides with u x , as can be easily checked. For the energies 
of solar neutrinos the matter effect on S13 is negligible. Furthermore, the oscillations of solar 
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neutrinos related to v 3 are averaged out. At the same time, v' e and v x form standard two 
neutrino system with vacuum mixing and mater effects. 

Using these features of dynamics, as well as projection of original flavor states onto 
z/, u x , v 3 one gets immediately the conversion probabilities: 


P* 

Pe. 

P,, 


c 13 Pee + s 13 
C 13^(l — Pee) 

Ci 3 (l — T 7 ) (1 — Pee) + S 13 C 1 3 (1 A Pee) i 


(37) 


where P ee is the //(-survival probability in the z/ — v x system: P ee = P ee ( 6 1 2 , A m\ 2 , V) and 
V = c\ 3 V e — V^(l — y). As in the previous case, the parameter r/ determines fraction of v s in 
the state v x to which v' e converts, and also it gives total fraction of the v s in the solar pair 
(9). 

Using (37), we get fluxes (in the units of the SSM flux) measurable by different reactions: 

&CC = fB [^13Pee + «l 3 ], (38) 

®nc = Jb[ 1 — (1 — P ee )r]cl 3 ], (39) 

= fB[c\ 3 Pee + 4 3 + rc w(l ~ Ae)(l ~ V) + rslii 1 + Pee)]- (40) 

Comparing with Eq. (13,14, 15) we find that corrections are indeed of the order s^ 3 . For 
the allowed values of s 33 these corrections have small impact on the data fit, as one can see 
from our analysis in [21]. 

The fluxes cc, ®es depend on two combinations of parameters, as in Eq. (17), 

with substitution x —> x\ y —> y', where 


x' = c* 3 x, y' = cl 3 + s 2 13 (x + f B ), 


(41) 


and x and y are defined in (16). So, the problem of degeneracy is not resolved with s 33 
corrections. Notice also that the fluxes in (38 - 40) satisfy the relation (18). 

The KamLAND signal is determined by the survival probability 

Pee = c\ 3 P2 + S 33 , (42) 

where P 2 is the two neutrino vacuum oscillation probability determined by A m\ 2 and # 12 - 
The factor cf 3 leads to additional suppression of the KamLAND signal and shift of the 
allowed regions to smaller # 12 (see detailed analysis in the 3// context in [34]). 


5.2 Small additional Am 2 

Let us consider the case when even small v e mixing beyond the solar pair can produce a 
strong effect. Suppose that, 
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(i) v e mixes with v x in the solar pair, v\ and v 2 , which has A m\ 2 and tan 2 d\ 2 in the 
LMA region, as in (35); 

(ii) the state zq ~ v' x (see (35)) is heavier than U\ and the mass splitting is A m 2 A1 ~ 10~ 6 
eV 2 . (The mass m\ can be very small or even zero). Furthermore, v e has small admixture 
in z/ 4 : tan 2 0 i 4 ~ 1CT 3 , so that v e — v' x oscillation parameters are in the SMA region. 

(iii) the admixture of v e in the z / 3 is zero (for simplicity). 

Then, similarly to (35) we can write the mixing as 


v\ = cosdx 2 (ci 4 V e -sx 4 v' x ) - sm6 12 v x , 

v 2 = sin 9 l2 (c 14 v e - s 14 z 4 ) + cos #12*4, 

F3 = c 13 v' a , 

v 4 = C141 4 + s 14 v e . (43) 


Here v' x is the combination of v a and v s orthogonal to v x ( 6 ), and v' a is the combination of 
v ^ and v r orthogonal to v a (7). 

The system has two resonances associated with A m\ 2 and Am 41 . Let us estimate 
qualitatively the conversion effect in different energy ranges. 

In the first resonance related to Am 2 2 the conversion is completely adiabatic. So, at 
the high energies, E > 7 — 8 MeV, where crossing of the resonance occurs the survival 
probability equals P^ e « sin 2 0 i 2 . At low energies, E < 0.5 MeV, (no 12-resonance crossing) 
effect driven by Am 2 2 is reduced to averaged vacuum oscillations. The result of conversion 
is similar to the one in usual three neutrino system: 


pL _ 
r ee ~ L 12 


Po + 


12 ; 


(44) 


where P 2 is the survival probability in the resonance related to A m\ A . It has a typical 
form of the survival probability of the SMA solution and characterized by Am 2 4 and mixing 
parameter tan 2 6 = tan 2 0i 4 / cos 2 Q\ 2 - Appearance of P 2 in this expression leads to an 
additional suppression of the pp-neutrino flux. 

In the intermediate energy range, E = (0.5 — 7) MeV, the conversion driven by Am 2 4 
will lead to flattening of the adiabatic edge of the suppression pit due to Am 2 2 . 


5.3 KamLAND result and LMA 

In the best fit point of the LMA region we predict Rkl = 0.66 for zero 77 [ 21 ], At 3 cr 
level Rkl < 0.73. In presence of the sterile component this bound can be slightly relaxed: 
Rkl < 0.76. For other solutions of the solar neutrino problem we predict Rkl ~ 1-0 if 
S 13 = 0. For non-zero S 13 the oscillations driven by A m 2 atm will lead the averaged oscillation 
result at KamLAND: Rkl = 1 — O. 5 sin 2 20 i 3 . For maximal allowed values of S 13 we get 
Rkl = 0.90 — 0.95 depending on A rn 2 atrn . What if Rkl = 0.76 — 0.90 will be found? 
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The following comments are in order: 

If solution of the solar neutrino problem is in the LMA region, a strong suppression is 
expected in KamLAND. In this case, introduction of additional sterile states will not help 
in the single Am 2 context, as we have found in the Sect. 3. 

Let us discuss what happens if two or more Am 2 contributes. 

Oscillation parameters can be taken beyond the LMA region in such a way that sup¬ 
pression in the KamLAND experiment is weaker. In this case, however, the conversion of 
solar neutrinos will not describe data in the two neutrino context. Then the parameters 
(mass, flavor mixing) of the z / 4 can be selected in such a way that conversion driven by z / 4 
will not change KamLAND result but correct description of the solar data. 

Such a possibility can be realized in the scheme described in the previous section 5.2 
with A m \ 2 ~ 6 ■ 10 -5 eV 2 . Indeed the KamLAND result is determined by the oscillations 
driven by Am 2 2 and 6 * 12 . For tan 2 6* 42 ~ 0.10 — 0.15 we get Rkl = 0.80 — 0.85. 

As far as the solar neutrino signal is concerned, the high energy signal will be suppressed 
by sin 2 9, so that large boron neutrino flux: J'b ~ 1.6 is required to explain the SNO and SK 
results. According to our previous consideration, in the single Am 2 context this point would 
be excluded by the Gallium experiment (large Qcie) and by SK (distortion of spectrum). 
The effect of Am 2 4 changes a situation: appearance of P 2 in expression for probability at low 
energies (44) leads to an additional suppression of the pp-neutrino flux, and consequently, 
Qce rate. In the intermediate energy range, E = (0.5 — 7) MeV, the conversion driven by 
Am 2 4 will relax the SuperKamiokande lower bound on mixing which follows from absence 
of spectrum distortion. The detailed analysis of this possibility will be given elsewhere [36]. 


Another possibility to reconcile an intermediate suppression of signal in KamLAND and 
explanation of the solar neutrino data is to split two problems. Suppose the solar pair has 
A m 2 2 in the LOW or VO regions, so that it does not produce any effect in KamLAND. The 
state V 4 which consists predominantly of the sterile component has the mass split: Am 44 > 
ICO 4 eV 2 . Then the suppression of the KamLAND signal is determined by admixture of v e 
in z/ 4 : Rkl = 0.5sin 2 2d 14 . For instance, Rkl = 0.85 can be achieved at sin 2 9 U = 0.09. 
Effect in solar neutrinos will be determined by P e s “ n = c\ A P ee + s 44 , where P ee is the survival 
probability for the conversion driven by A m\ 2 . Additional factor c 44 ~ 0.8 can even improve 
the fit of the data for LOW and SMA [21]. 


6 Conclusions 

1. We have studied properties of mixing of the sterile and active neutrinos. In the single 
Am 2 context (when electron neutrino is present in two mass eigenstates only) the problem 
is reduced precisely to two neutrino problem with mixing of v e with v x in the solar pair (y x 
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is the combination of the active and sterile neutrinos). 

2 . We have performed the global fit of the solar neutrino data in the single Am 2 context, 
treating the boron neutrino flux as a free parameter. The best fit corresponds to zero 
admixture of sterile component, r) = 0. We get the upper bound r/ < 0.26 (0.64) at la (3a). 
With increase of q the best fit point shifts to smaller Am 2 and tan 2 0 i 2 . 

3. Due to degeneracy of parameters, no one single observable gives the bound on rj. The 
bound on sterility appears as an interplay of the SNO results whose allowed region of the 
oscillation parameters shifts to smaller mixings and Am 2 with increase of rj and Gallium as 
well as SK results which restrict the mixing and Am 2 from below. 

4. Further more precise measurements of the NC event rate at SNO will strengthen the 
bound: rj < 0.5 (3cr). 

5. Implications of KamLAND depend on specific results of KamLAND measurements. If 
KamLAND result (rate, spectrum) corresponds to the best fit point of the solar neutrino 
analysis, then after 3 years of the KamLAND operation the limit will be strengthened down 
to r) < 0.19 (0.55) at la (3cr). The limit will be stronger if the KamLAND best fit point 
will be at larger tan 2 d 12 and Am 2 than the solar neutrino data give at present. In contrast, 
if KamLAND selects smaller tan 2 0 12 and Am 2 , the limit will be weaker, and moreover, the 
data may prefer non-zero sterile admixture. 

6 . The presence of sterile neutrinos (in the single Am 2 context) relaxes the upper bound 
on the predicted rate Rkl only weakly from 0.73 to 0.76 (3cr). Larger values of Rkl can be 
reached in the two Am 2 context. 

7. If electron neutrino mixes also in other mass eigenstates, the effects depend substantially 
on values of other Am 2 . 

If additional A m\ A 3> A m\ 2 the small corrections appear to formulas with a single Am 2 
which are proportional to the mixing of u e in additional mass eigenstates, e.g., s 2 3 . If Am 2 , 
is in the MSW range an additional mixing can be resonantly enhanced in matter producing 
much larger effect. In this case intermediate values of Rkl can be obtained. 
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tan 2 0 


Figure 1: The best fit points and 3a allowed regions of oscillation parameters for different 
values of rj: rj = 0.0 (solid line), 0.2 (long-dashed line), 0.4 (dashed line) and 0.6 (dotted 
line). Boron neutrino flux is taken as free parameters. Also shown lines of constant Jb 
which minimize % 2 . 
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Figure 2: The allowed regions in fs — V plane, at la, 90%, 95%, 99% and 3a C.L.. The 
best fit point (marked by star) is at q — 0. 
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Figure 4: The la allowed region (right upper parts) from the analysis of the Su- 
perKamiokande data for different values of 77 . Also shown are lines of constant fs which 
minimize y 2 fit of the SK data for different 7 ]. 
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Figure 5: The predicted value of <f> 7 vc as a function of rj. The band corresponds to the lex 
region. The horizontal lines show the present SNO result. 
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0.6 0.65 0.7 0.75 0.8 


Figure 6: a) Minimum \ 2 obtained from analysis of the solar neutrino data as a function 
of the KamLAND rate, Rkl > V — 0 (dashed line) and for rj taken as a free parameter 
(solid line), b) Value of r/ (solid line and left axis) and Jb (dashed line and right axis) that 
minimize the x 2 presented in the panel a). 






True values for the neutrino parameters that are indicated by a stars, b) The dependence 
of the x 2 on V f° r these points from a combined analysis of the simulated KamLAND and 
the solar neutrino results. 
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tan 2 0 

Figure 8: The upper panel: the lines of constant lcr upper bound (figures at the curves) on 
r] which can be obtained from from combined analysis of simulated KamLAND and solar 
neutrino results. The lower panel: the same as in the upper panel but for la lower bound 
on r). 













